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ABSTRACT 


From the linearized kinetic theory of stimulated cyclotron 
radiation, and from related physical arguments, some design cri- 
teria are deduced for spiraling-beam cyclotron-resonance oscilla- 
tors for possible use at millimeter-wave frequencies. While an 
optimum choice of the geometry of the interaction region is 
important, the most essential design problem is the control and 
minimization of the axial velocity spread on the . spiraling beam. 
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SUMMARY 


Recent work on periodic-beam microwave devices has indicated 
that these devices show some promise for development into power- 
ful, tunable sources in the millimeter and submillimeter wave- 
length range. However, no comparative evaluation of different 
designs has been attempted to date. In this report, some design 
criteria are deduced for spiraling-beam cyclotron-resonance 
oscillators, based on the linearized kinetic theory of stimulated 
cyclotron radiation and on related physical arguments. In the 
absence of limiting-amplitude calculations, the start-oscillation 
current is used as an indicator of device activity and, hence, a 
qualitative predictor of efficiency. It is shown that while an 
optimum choice of the geometry of the interaction region is impor- 
tant, the most essential design problem is the control and mini- 
mization of the axial velocity spread of the spiraling beam. 


INTRODUCTION 

Cyclotron resonance oscillators are devices where an ensemble 
of mono-energetic electrons, orbiting in a uniform magnetic field, 
are stimulated by an rf electric field to emit coherent radiation. 
This is not to be confused with the spontaneous (and generally 
incoherent) radiation of accelerated classical charges (see Appen- 
dix A) . Several different devices of this class have been demon- 
strated recently, both with spiraling electron beams (refs. 1-5) 
and with cycloidal sheet beams (ref. 6) . These devices show some 
promise for development into powerful, tunable sources in the 
millimeter and submillimeter wavelength range. As such, they 
would help to open up a range of the electromagnetic spectrum 
that until now has not been very useful for communications because 
of a lack of powerful sources. 

The purpose of this report is to serve as a summary of work 
done to date on an in-house project at the Electronics Research 
Center, which is aimed at a broad engineering evaluation of cyclo- 
tron resonance devices. In particular, it is desired to assess 
the inherent limitations, if any, on power output, efficiency, 
and frequency, and in the process to derive information that will 
result in optimum design. The theoretical work is complemented by 
experimental data designed to verify certain predictions and im- 
prove the understanding of the physical phenomena taking place. 
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axial beam velocity in the interaction region 

transverse beam velocity in the interaction region 

corkscrew pitch 

cyclotron frequency 

charge-to-mass ratio 

permeability of free space 

current 

corkscrew radius 

axial beam velocity at some point along the axis 

transverse beam velocity 

linear beam diameter at cathode 

linear beam diameter after compression 

axial magnetic field linking the cathode 

axial magnetic field at some point along the axis 

transverse magnetic field 

total beam diameter 

cyclotron radius 

ac current density 

rf electric field vector 

axial and transverse components of the wave vector 
dc beam current 

radian frequency of rf cavity fields 
energy stored in cavity 
loaded Q 
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transit time 
unloaded (internal) Q 
external Q 

equivalent shunt resistance of cavity 
equivalent shunt conductance of cavity 
equivalent electronic beam-loading conductance 
rf voltage drop along electron trajectory 
beam clearance in cavity 
relativistic cyclotron frequency 
electron rest mass 

relativistic mass change parameter 


REVIEW OF PREVIOUS WORK 

Cyclotron resonance oscillators have had a variety of 
experimental embodiments, and some of these designs are sum- 
marized briefly in Table I. As might be expected, the design of 
a prototype model in connection with the in-house work, drew 
heavily on the information contained in Table I. Two design 
decisions are basic: the choice of the beam trajectory, and the 

choice between resonant or backward-wave interaction. A simple 
helical beam, generated by a magnetic corkscrew, was chosen for 
its relative simplicity and precise controllability in the 
laboratory . 

From Table I, it is evident that devices based on backward- 
wave interaction have generally shown greater efficiency than 
resonant-cavity devices. However, the former require a longer 
region of homogeneous magnetic field. Thus, in the absence of 
any analytical evidence that backward-wave devices are inherently 
more efficient, it was decided for the in-house experiment to 
take advantage of the economy in magnetic-field volume offered 
by the resonant-cavity approach. 

The basic operation of a helical-beam cyclotron resonance 
oscillator is described in Figure 1. An electron beam is formed 
by a triode gun that is immersed in an axial magnetic field of 
two to three times the Brillouin value. A small transverse 
magnetic field component forces the electrons into slightly 
spiraling trajectories, and the beam then drifts axially into a 
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TABLE I 

COMPARISON OF EXPERIMENTAL PARAMETERS FOR SOME CYCLOTRON RESONANCE OSCILLATORS 


Experiment 

Bott 
(ref. 1) 

Hirshf ield 
Bernstein and 
Wachtel 

(ref. 2) 

Gaponov , 
et al . 

(ref. 4) 





Type of gun 

Pierce gun 
immersed 

Parallel -plane 
diode, immersed 

Conical 

cathode 

Electron 

trajectory 

Spiral 

i 

Spiral 

Spiral 

Beam voltage 

20 kV 

5 kV 

19 kV 


Schriever and Antakov, 

Johnson et al . 


Johnson 
(ref. 3) 



Cathode current 

Type of 
trajectory 
modulation 


320 mA 


Spiral 


10 kV 


1300 mA 


(ref. 6) 


Kino gun 


Cycloid 


14 kV 


600 mA 


■ 

I 


Mirror ratio 


Type of 
interaction 


Power 

Frequency 

Efficiency 1 ' 


Corkscrew 


Combination 

resonant and 

backward-wave 

Circular TE 

mn 

1 W 

143 GHz 


Corkscrew 


Resonant , 
circular 


Crossed-f ield Crossed-f ield Crossed-f ield 

injection injection injection 


Resonant , 
circular 


TE 011 

TE 021 


10 mW 

190 W 

720 W 

) . 8 GHz 

25 GHz 

(2nd harmonic) 

9 GHz 

_ 

3% 

6% 



J 


Backward-wave Backward-wave 


850 W 
38 GHz 


Kulke and 
Wilmarth 

(ref. 7) 


Pierce gun 
immersed 


Spiral 


15 kV 


25 mA 
Corkscrew 


Resonant , 
rectanaular 
TE. 1 . 
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region of increasing magnetic field, where more and more axial 
kinetic energy is transformed into transverse energy. Once the 
beam reaches the point where nearly all its momentum is in the 
transverse direction, it drifts into a microwave cavity that is 
resonant in a TE mode slightly above the cyclotron frequency of 
the orbiting electrons. Within a narrow range of drift velocity 
and magnetic field, this resonance will couple to the spiraling 
beam, and gain or oscillations are obtained. The spent beam is 
then collected. Power is coupled from the cavity through an iris 
into the waveguide output. 

Throughout this report, the small-signal interaction between 
the spiraling beam and the electromagnetic field is discussed in 
terms of the linearized kinetic theory used by Hirshfield, 
Bernstein, and Wachtel (ref. 2) . Other approaches are possible 
and have been pursued elsewhere. One of these is the Pierce- 
type analysis (ref. 8) which was used, for example, by Chow and 
Pantell (ref. 9) in their analysis of a traveling-wave spiraling- 
beam device. In this method, one describes first the rf modula- 
tion on the beam due to the circuit fields, then the circuit 
fields induced by the beam perturbation, and finally one requires 
the two expressions to be consistent. The resulting dispersion 
relation is formidable, but with the aid of a high-speed computer 
it can be solved numerically for the complex propagation constants. 
The small-signal gain and the start-oscillation conditions can 
then be calculated by imposing appropriate boundary conditions. 

Two different physical mechanisms have been identified which may 
cause the necessary bunching (phase-focusing) of the orbiting 
electrons into a decelerating phase in the transverse electric 
fields of the rf cavity. Chow and Pantell based their analysis 
on one such mechanism where the dc rotation of the electrons and 
the transverse ac magnetic field of a traveling wave produce a 
z-directed Lorentz force, and this tends to bunch electrons into 
the decelerating phase of the traveling transverse electric field 
vector of the waveguide. This process will hereafter be called 
"Pantell bunching." However, a phase-focusing effect can also 
arise from the relativistic mass change of the orbiting electrons. 
If the orbiting electrons are initially assumed to have random 
phases, then those with initial phases such as to gain energy 
from the rf fields become heavier and accumulate phase lag, and 
those with initial phases such as to lose energy to the rf fields 
become lighter and accumulate phase lead. Sehn (ref. 10) re- 
derived the dispersion relation applicable to backward-wave 
interaction in relativistic form, and he concluded that the 
relativistic bunching mechanism is roughly as effective as the 
Pantell bunching, and, therefore, cannot be neglected. 

In recent experiments at the University of Utah (ref. 11) , 
both forward and backward traveling-wave interaction have been 
observed, and were clearly distinguishable in terms of frequencies. 
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In addition, a third mode of interaction occurred at a frequency 
halfway between the other two modes, which was labeled the 
cyclotron resonance mode, but this mode was not analyzed further. 

Although the Utah group did not consider any relativistic 
bunching effects, their theory did allow them to identify accu- 
rately beam interaction with the forward or backward traveling 
waveguide modes. Their theoretical prediction of start oscilla- 
tion current has not been verified accurately by experiment, 
however, and thus it is difficult to tell at this time what 
quantitative effect, if any, their neglect of relativistic effects 
does have. A clear confirmation or refutation of Selin' s conclu- 
sion, therefore, does not yet exist. 

It is of interest to compare the theoretical synchronism 
conditions that apply to traveling-wave and standing-wave inter- 
action. Interaction between the beam spiraling with gyro- 
frequency, ft, and the forward-traveling circuit wave of fre- 
quency, go, requires that 


a) = ft + k || u , 


( 1 ) 


where k|| and u are the axial propagation constant and the axial 
beam velocity, respectively (ref. 11). As will be shown later in 
the section discussing the gain function, an almost identical 
synchronism condition exists for the standing-wave interaction, 


a) = ft - xk || u ( 2 ) 

where x - (-0.85) corresponds to the desired condition of peak 

negative beam loading of the cavity. The standing-wave inter- 
action frequency, therefore, would be expected to differ only 
slightly from that for the forward-traveling wave interaction. 

One is led to the intuitively appealing conclusion that the 
standing-wave interaction actually takes place selectively with 
one of the component traveling waves. An experimental confirma- 
tion of this conclusion, however, is not yet available. Tied in 
with this problem is the unanswered question of whether traveling- 
or standing-wave interaction leads to a basically more efficient 
device . 

In Pierce's original analysis of beam-wave interaction 
(ref. 8), the coupling between the normal modes of a beam and of 
a waveguide is considered implicitly, as certain propagation 
constants are dropped from the dispersion relation. A more 
explicit coupled-mode formalism was developed later, first for 
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a linear beam (ref. 12) and then for a beam executing transverse 
motion under the influence of rf forces (ref. 13) . 

This coupled-mode formalism has been applied extensively, 
for example, by the Utah group (ref. 11) in their thorough 
follow-up to Chow and Pantell's work on transverse traveling wave 
devices . 

Mclsaac (ref. 14) has derived a relativistically correct set 
of coupled-mode equations, based on the six normal modes of a 
beam with spiraling dc motion. This work also is still in prog- 
ress, and no correlation with experimental data has so far been 
attempted. 

An important advantage of the coupled-mode formalism is its 
capability to identify possible beam-circuit interaction modes. 

To proceed from the coupled-mode expressions to tractable engi- 
neering formulas, e.g. , for the start-oscillation current, 
however, seems to be relatively difficult. 

A somewhat different approach was taken by Hsu and Robson 
(ref. 15) who analyzed the spiraling-beam interaction with a 
Cuccia coupler by assuming a small rf perturbation on the dc 
electron velocity vector and carrying this directly through the 
non-linear, relativistic equation of motion. This approach was 
later extended by Hsu (ref. 16) to include beam interaction with 
the transverse electric fields of a resonant cavity. With suit- 
able assumptions, Hsu was able to get good agreement with the 
start-oscillation currents measured earlier by Bott (ref. 1). 
However, these assumptions involve the ratio of transverse to 
axial velocity which, although crucial to the result, is very 
difficult to verify. experimentally. 

In fact, the lack of accurate knowledge of the axial veloc- 
ity and of the velocity spread in the interaction region appears 
to plague all predictions of start-oscillation current, regard- 
less of the theory on which they are based. In this sense, the 
distinctions between the different theories become almost irrele- 
vant, because any of them only provide rough predictions of the 
start-oscillation current. 

Further experimental work which should also be mentioned is 
the device built by H. Jory (ref. 17), where an initially linear, 
0.100 inch diameter, beam is injected into an X-band accelerator 
cavity driven with up to 75 kW of pulsed power . The beam emerges 
from this cavity in a spiraling trajectory with a kinetic energy 
of 300-500 keV (mainly transverse) , and then enters axially into 
a cylindrical TE n ,1,1 ("barrel") resonator. As the beam spirals 
through the cavity, almost grazing the wall, it will interact 
with the tangential component of the periodic rf electric field 
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patterns in what closely approximates conventional linear-beam 
periodic-field type interaction. The combination of the very 
high beam velocity and the overmoding of the cavity (up to n = 14) 
has enabled the device to operate at nearly 10 times the electron 
cyclotron frequency, with outputs up to 6 W at 84 GHz. The 
spiraling trajectory of Jory's beam could lead one to attempt an 
analysis of his device in terms of Pantell-type or relativistic 
bunching, but the approach followed by Jory, in terms of con- 
ventional bunching along the circumference, is clearly more 
appropriate. Thus, this tube really is not a cyclotron resonance 
device . 

As has been mentioned previously, the small-signal kinetic- 
theory approach as used by Hirshfield, Bernstein, and Wachtel 
(ref. 2) will be followed in this report. This approach is 
borrowed from the realm of plasma physics, as it involves solving 
the collisionless Boltzmann equation, and it is definitely 
limited to the small-signal regime, e.g., the prediction of 
electronic gain and of start-oscillation currents. It is attrac- 
tive because of its algebraic compactness and because it auto- 
matically takes into account both the relativistic and the 
Pantell-type bunching. Also, the kinetic-theory approach in 
principle yields solutions for any given electron velocity dis- 
tribution. This latter capability is important, as both theory 
and experiment indicate that the width of this distribution con- 
stitutes perhaps the most important design characteristic of a 
cyclotron resonance oscillator. 

Some problems encountered in generating a beam having the 
desired narrow velocity distribution are discussed in the follow- 
ing section. 


FORMATION OF A SPIRALING ELECTRON BEAM 

The choice of beam geometry and the manner in which the 
beam is generated have a profound effect on the operation of 
cyclotron resonance devices because both affect the velocity 
distribution of the electrons interacting with the cavity fields. 
In Table I, a comparison has already been made between the 
pertinent design data for some experimental devices. A tremen- 
dous amount of work has been done by many authors in the field of 
electron optics, and no attempt will be made here to provide a 
comprehensive evaluation of different gun systems suitable for 
spiraling-beam oscillators at a given power level and frequency. 

The design to be discussed here already has been mentioned 
in connection with Figure 1. The transverse perturbation is 
of the order of one percent of the axial field and is set up by 
a "corkscrew" magnetic field (ref. 18) that is generated by a 
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current-carrying bifilar helix wound externally around the non- 
magnetic, evacuated drift tube. The length of this helix is 
somewhat arbitrary, but there is an obvious trade-off between 
the length and the helix current. The transverse magnetic field 
vector is screw-symmetric along the beam and thus forces the 
electrons into slightly spiraling trajectories with pitch identi- 
cal to that of the corkscrew, and the fraction of axial kinetic 
energy so converted into transverse energy will depend on the 
corkscrew length and on the current. The corkscrew pitch p c 
follows from the synchronism requirement 


Pc (z) _ 2n 
u(z) io b 


(3) 


where u(z) = axial beam velocity, a function of distance, and 
w b = ( e / m ) B l = cyclotron frequency. As the axial beam velocity 
in many cases is only a weak function of distance, the corkscrew 
may be wound with constant pitch rather than with the tapered 
pitch that is called for by Eq. (3). For example, a total change 
in axial kinetic energy of about 10 percent corresponds to a 
change of only 5 percent in velocity over the total length of the 
corkscrew. For a 3-4 turn corkscrew, the required tapering would 
call for a winding accuracy beyond the mechanical tolerance 
limits within which each turn can be placed. However, one must 
realize that a nonideal corkscrew (i.e., one with fringing end 
fields and with imprecisely spaced turns) is bound to degrade 
even an incoming mono-energetic beam so that some velocity spread 
will appear at the corkscrew exit. Similarly, a moderate spread 
in the velocity of the incoming beam is, in effect, amplified 
during passage through the corkscrew region. Analytical expres- 
sions for the geometry of an optimally efficient corkscrew, and 
the current required to produce a given energy conversion with a 
given beam, are available (ref. 19) . The required calculations 
are laborious , but a computer program has been written and the 
results have been verified experimentally with great accuracy.* 

As the corkscrew current usually can be adjusted over some range, 
the following rough approximation may be used if the desired 
transverse field is known. The bifilar helix is modeled by a 
series of interleaved planar current loops (Figure 2) . The axial 
field components cancel, and the transverse field becomes that of 
two straight parallel wires spaced a distance 2 r c apart. On the 
beam axis, one then finds for the transverse field strength 


B 1 “ V Ar c 


(4) 


*W. Getty, University of Michigan, to be published. 
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Figure 2.- Approximation used 
for calculating the trans- 
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For example, on a corkscrew with a 2-cm mean diameter and a 
10-cm pitch, Eq. (4) was found to predict the measured field 
perturbation (0.43 Gauss per Ampere) within better than ten per- 
cent . 


The radial decay of the transverse magnetic field perturba- 
tion away from the corkscrew has a modified Bessel function 
dependence (ref. 18). 



In order to make f (r) - constant over the beam cross-section, 
it is clearly desirable to make the corkscrew pitch p c and hence 
also, the beam velocity u large. 

One might note here that immersed flow implies essentially a 
slipping-stream beam, in which the axial velocities are non- 
uniform over the beam cross section, in contrast to a Brillouin 
beam. Consequently, the synchronism condition Eq. (3) cannot be 
met by all electrons simultaneously, so that the inherent velocity 
spread on the incoming beam will, in fact, be amplified by the 
corkscrew action. To avoid this, one clearly should either use 
Brillouin flow or else resort to a very thin beam. 

In addition to the magnetic corkscrew, other methods of 
shaping the initial beam trajectory are possible, such as 
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electrostatic corkscrews or deflection plates, and various ways 
of injecting the electrons at an angle to the magnetic field 
lines. Strictly as a laboratory device, the simple magnetic 
corkscrew seems to perform rather well, however. 

Once started on their spiraling trajectory, the electrons 
are allowed to drift into a region of increasing magnetic field 
Bj_, here called a mirror region. It is assumed that the change 
in field is slow enough so that the electrons will not cross flux 
lines, and the magnetic moment will thus be conserved (the 
adiabatic assumption). Then it can be shown (ref. 20) that, as 
increases, the orbits decrease in size and also move closer 
together, each guiding center remaining on a flux line and each 
orbit within a tube of flux. The kinetic energy of the circular 
motion will increase directly proportional to the magnetic field 
strength, at the expense of the axial energy (see Figure 3) , 


magnetic moment 


A 2 mW 



constant 


(6) 


MAGNETIC FLUX LINES 



Figure 3.- Magnetic compression of a beam of spiraling 
electrons, assuming the compression to be 
adiabatic. 


A typical beam cross-section such as might result with a 
15-keV beam after nearly 100 percent conversion of axial energy, 
orbiting in a 3500-G magnetic field, is shown in Figure 4. 
Throughout this discussion, the assumption is being made that 
the electrons are spiraling individually, each about its own 
guiding center. While this assumption is certainly true of a 
linear, confined-flow beam, direct experimental observation of 
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A) ELECTRONS SPIRALING INDIVIDUALLY 



B) ELECTRONS SPIRALING COLLECTIVELY 



Figure 4.- Cross-section of 15 keV beam of spiraling electrons 
after nearly 100 percent energy conversion. 


spiraling beams has shown that under the influence of the cork- 
screw perturbation, the beam as a whole will assume a helical 
shape, fixed in space. However, after passage into the mirror 
field, with nearly 100 percent of the axial kinetic energy con- 
verted, the pitch of the helical electron orbits will approach 
zero in either case, and the outer diameter of the resulting 
rotating electron cloud will ideally be the sum of the original 
beam diameter and the orbit diameter. An experiment was carried 
out in which the beam was allowed to impinge on a tungsten-mesh 
target, exciting the latter to local incandescence. The incan- 
descent outline of the beam was viewed through a glass window 
while the pitch of the electron orbits was gradually decreased 
by increasing the corkscrew current. The fact that the beam 
outline then was seen to bloom out symmetrically from the center, 
without any observable lateral movement across the target, 
clearly favors the picture of electrons spiraling individually 
at this point. 

Actually, with the zero-space-charge conditions assumed in 
the small signal kinetic theory below, the distinction between 
individual and collective spiraling is unimportant, because the 
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electrons then act in any case as independent, identically phased 
sources driving the cavity (ref. 7). In this report, for sim- 
plicity, individual spiraling of the electrons will be assumed. 

With this assumption, one can calculate the beam diameter 
as a function of the axial magnetic field strength. The initial 
linear beam diameter d^ will undergo magnetic compression. 


d 


1 



(7) 


where Bq = magnetic field linking the cathode, and d^ = linear 
beam diameter at some distance along the axis, with B = B^. 

Using the data in Figure 4, one superimposes the cyclotron orbits 
of the individual electrons on the compressed linear diameter, 
in such a way that each electron retains its original guiding 
center. The resultant overall beam diameter then becomes 


D = d^ + 2r £ 


(8) 


where rj-, = cyclotron radius. A plot of D versus axial distance 
is shown in Figure 5, for the device described in Figure 1. A 
beam of 15-keV electrons that initially has a 2-mm diameter, 
first expands while passing through the corkscrew and then con- 
tracts under the influence of the increasing axial magnetic 
field, reaching a final diameter of 3 mm at B = 3.43 kG (cyclo- 
tron frequency 9.6 GHz), with all its kinetic energy in the 
transverse direction. Note that the beam diameter shown in the 
corkscrew region is only approximate. 

We now return to the problem of velocity spread. As 
mentioned before, winding irregularities and the transverse field 
variation on the corkscrew are apt to generate some velocity 
spread on an initially mono-energetic beam, or will increase 
whatever velocity spread already exists on the beam. In addition, 
the action of the mirror field will itself be such as to increase 
the relative width of any velocity spread. The conversion of 
axial to transverse energy in effect moves the initial distribu- 
tion which may have been, for example, a narrow Gaussian centered 
about 15 keV, closer to the origin, and the distribution may be 
truncated or otherwise distorted as the slowest electrons begin 
to be turned around by the magnetic mirror. Also, space charge 
effects will begin to be important as the beam is slowed axially, 
and this will result in potential depression near the beam center, 
and will thus increase the velocity spread. 
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une should note that these effects all arise even with a 
mirror field that completely satisfies the adiabatic assumption, 
i.e., where the spatial rate of change of magnetic field is very 
gradual. "Gradual" has been described by Allis (ref. 20) with 
the inequality 



(9) 


which means, essentially, that any spatial change in should 
be distributed over a distance large compared to the cyclotron 
radius. For example, at the steepest point of the field profile 
shown in Figure 1, the characteristic distance is 


B 

grad B 


90 mm , 


and this is far greater than the cyclotron radius at that point, 
r^ = 1.64 mm. 

If the adiabatic assumption is not satisfied, then the 
electrons will begin to cross flux lines, and this will further 
increase the velocity spread. As an air-core solenoid system 
will generally result in smaller field gradients than one where 
iron pole pieces are used, the former is preferable and thus has 
been used with the in-house X-band device. For devices that 
employ superconducting magnets, pole pieces would saturate and 
thus are not used in any case. 

A word should be added here about standing-wave (resonant) 
versus traveling-wave (backward-wave) interaction with the 
spiraling beam. It has been mentioned that the resonant inter- 
action is more economical in terms of the required volume of 
uniform magnetic field. However, the existence of velocity spread 
on the beam probably favors axially long interaction regions, 
because these, together with some inevitable magnetic field non- 
uniformities, may well serve to utilize a greater velocity range 
of electrons for the interaction. This may in fact be the reason 
why the devices using backward-wave interaction have generally 
been found to be more efficient than standing-wave devices. More 
definite information should become available once a planned 
numerical large-signal analysis has been completed. 

One more question which must here be left unanswered, is 
how much of the conversion of axial to transverse energy should 
be accomplished by the corkscrew, and how much by the magnetic 
mirror. Resolution of this problem will become clearer once 
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more is known about the way in which corkscrew and mirror fields 
operate on the velocity distribution. 


START-OSCILLATION CURRENT FOR A MONO-ENERGETIC BEAM 

In this section, the small-signal theory will be summarized. 
The electron beam inside the cavity is considered as a cloud of 
slowly drifting electrons, each orbiting at the cyclotron fre- 
quency about a magnetic flux line. The inclusion in the theory 
of an arbitrary velocity distribution in principle allows one 
to take into account varying degrees of velocity spread. The 
electron cloud is assumed to be made macroscopically neutral by 
the presence of immobile ions, and is assumed to be sufficiently 
tenuous such that electron-ion collisions and electron-electron 
collisions (rf space-charge effects) can also be neglected. Given 
the original electron velocity distribution, the linearized 
Boltzmann-Vlasov equation is solved for the perturbed (i.e., ac) 
part of the distribution function in terms of the known electro- 
magnetic field components within the cavity. The ac density 
follows by integrating over the electron velocity distribution, 
and the power flow from the electrons to the rf cavity fields 
is derived by integrating the product over the cavity volume. 

(It is assumed that the cavity modes are only insignificantly 
perturbed by the presence of the free charges.) The result is, 
after some manipulation, 


p = - j - E 2 ~ P (kj_a) jdwjdu f Q (u,w) ^ G (jj (x)[l+ 6Q (x) ] 

mk I, J J 

( 10 ) 

This equation describes the following situation. The column of 
electrons, with radius "a," drifts axially through the cavity, 
with the rf electric field polarized at right angles to the beam. 
The electrons initially have a velocity distribution fg(u,w), 
where u and w are the axial and transverse velocity components, 
respectively. The term p (kj_a) represents the number of electrons 
per unit length of the beam, weighted by the variation of the rf 
electric field over the beam cross-section, and G w (x) [ 1 + 3Q W (x) ] 
is the so-called gain function, which describes the dependence 
of the absorption behavior on the transit time and on the synchro- 
nism conditions of the beam electrons. The remaining symbols 
have their usual meaning, 

e,m = electronic charge, mass 

kj| ,kj_ = axial and transverse components of the wave vector k 

E = amplitude of the rf electric field vector inter- 

acting with the beam. 
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In order to derive a simple approximate expression for the start- 
oscillation current, the following restrictions are now placed on 
Eq. (10). First, we assume a mono-energetic beam. 


f 0 (u,w) 



6 (u - Ug ) 6 (w - Wg ) 


( 11 ) 


The second assumption concerns p (k^a) • For the referenced ex- 
periment (ref. 7), a 3-mm diameter, spiraling beam traverses a 
rectangular TEjq^ cavity (f = 9.463 GHz), of length L = 2.14 cm 
and cross-section 2.29x1.02 cm, at right angles to the electric 
field vector, and the electric field variation over the beam 
cross section is very small. Hence, one may approximate p (kj_a) 
by a constant 


p 0 


u 0 e 


( 12 ) 


where I n = dc beam current. Substituting Eqs . (11) and (12) into 

( 10 ), 




P 




G W (1 + 6Q «> ) 


(13) 


For oscillations to start, the power extracted from the beam is 
just balanced by the losses in the cavity walls plus the external 
load. 


P = “2 

loss Q l 


(14) 


where U = stored energy and Q L has its usual meaning. Finally, 
equating Eqs. (13) and (14) , one finds the start-oscillation 
current 


I start 



o)U 1 

QLE 2*G aJ (l+ 3Q W ) 


(15) 
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where r = L/uq denotes the transit time through the cavity, and 
ky = tt/L was used, corresponding to the TE^oi resonance in a 
cavity of length L = 2.14 cm. 

As will be shown, Eq. (15) leads to results that are in 
reasonably close agreement with experimental data. Its real 
significance, however, is that it points out those physical 
parameters which bear on the small-signal electronic gain of the 
device and which probably also have some influence on the satu- 
rated power output and the efficiency. These factors will now 
be considered in greater detail. 


CAVITY LOADING AND GEOMETRY 

The inverse dependence of Istart on Ql was well confirmed 
through measurements (ref. 7) where Ql was varied by changing the 
external loading, with 




1 


Q 


ext 


(16) 


For convenience, these results are reproduced in Figure 6. 

It is illuminating to consider a simple model where the 
presence of the beam causes a negative conductance G e to appear 
in parallel with the cavity shunt conductance G s ^, and the 
start-oscillation condition corresponds to | G e | >. G s h» As one 
can write 





(17) 


where the first two terms depend only on the geometry of the 
cavity and the coupling iris respectively, it is clear that 


G 


sh 



and hence , 


(18) 




beam 


(19) 
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Figure 6.- Start-oscillation current versus 1 /Ql, 
for the device shown in Figure 1. 


at start-oscillation. One concludes that in the small-signal 
regime, the negative beam- loading conductance G e is directly 
proportional to the beam current, that is, the orbiting electrons 
act as a linear superposition of individual, identically phased 
sources driving the cavity fields. This agrees with the zero- 
space-charge assumption made earlier. 


In Eq. (15), the influence on I s tart °f the cav ity geometry 
is described by the factor e2/ojU. For want of a better term, 
this factor will be called the "shape parameter" of the cavity, 
because it provides a measure of the effectiveness of a given 
amount of rf stored energy, U, in setting up the interaction 
field intensity, E. Clearly, if a small Istart is desired, the 
cavity should be designed in such a way that most of the electric 
energy is stored in the region traversed by the beam. 


The similarity of 
shunt impedance of the 


the "shape parameter" 
cavity is evident. 


to the normalized 



( 20 ) 


i 
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where V stands for the peak value of the rf voltage drop experi- 
enced by an electron in transit through the cavity; i.e., the 
line integral of the rf electric field vector along the path of 
the electron. It is patently difficult to write such an integral 
for an electron orbiting in the electric field, especially since 
the number of orbits is not easily determined. However, for many 
practical cavities one may assume a reasonably constant electric 
field in the interaction space, and then arbitrarily let 


V = E 0 b (21) 


where b = gap spacing, determined mainly by the required beam 
clearance. With this assumption one then finds 


E 2 _ 2 R sh 

wU b 2 Q 0 


( 22 ) 


An approximate calculation of the beam diameter has already 
been carried out and is shown in Figure 5. 

Inside the cavity, the beam diameter is expected to be about 
3 mm, however, a larger gap spacing (7.62 mm) was assumed arbi- 
trarily to make allowance for some error in beam alignment. 

Figure 7 shows the calculated values of the "shape parameter" for 
some simple cavity geometries. All the applicable cavity dimen- 
sions are optimized to give maximum E^/wU at fo = 9.6 GHz, with 
the beam clearance as stated above. 

It is seen that the cylindrical TMqio anc ^ t ^ ie rectangular 
TEioi cavities have almost identical shape parameters. However, 
the latter has slightly smaller internal losses (greater Qo) and 
also is fitted more easily to a rectangular feeder waveguide. 
Hence, for operation at 9.6 GHz, the rectangular cavity is the 
best of the simple shapes; further improvement may be possible 
by going to more complex, empirical designs, e.g., by adding a 
center post to the cavity. 

A similar calculation may be made at 94 GHz. Assuming the 
same initial beam diameter, the final diameter will be 0.442 mm 
in an axial field of 33.6 kG, corresponding to a 94-GHz cyclotron 
frequency. This beam fits easily through a rectangular TEioi 
cavity consisting of a resonated section of WR-10 guide, with 
cross-section of 2.54x1.27 mm. 

At higher frequencies, the internal cavity losses will dis- 
sipate an increasing proportion of the total rf power extracted 
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Field Configuration 


Description 


E m/wU 



Rectangular 


(Square Base) 


Cylindrical TMq jq 


Cylindrical TE.-^ 

(Min . beam clearance 
does not apply) 


Cylindrical TEq^ 


(Min. beam clearance 
does not apply) 


fS 9 

4.04x10 ohms/m 


Q = 8400 
o 


4.08x10^ ohms/m^ 


Q = 6900 
x o 


6 9 

1.56x10 ohms /m 


0.34x10^ ohms/m^ 



Figure 7.- Shape parameter for some simple cavity geometries, 
optimized at f = 9.6 GHz. The beam path is indica- 
. ted by the dotted arrows; the electric field lines 
by the solid arrows. 
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from the beam. For example, the theoretical Qo for a rectangu- 
lar TE^oi cavity is 7800 at 9.6 GHz, and drops to 2800 at 94 GHz, 
varying roughly as (f)"i/2. The choice of a larger volume-to- 
surface ratio will increase Qq , as for example a spherical 
biconical cavity, heavily overmoded, can easily produce a Qq 
greater than 10,000 at X-band (ref. 21). However, Qo is then 
traded off against the shape parameter because in the overmoded 
cavity most of the electromagnetic energy will be stored away 
from the interaction region. The same comment can be made about 
quasi-optical (overmoded) resonators in general, including Fabry- 
Perot etalons. 

On the other hand, one application where small cavities have 
been used successfully even up to 700 GHz is that of frequency- 
multiplier klystrons (ref. 22), and this is one example which 
illustrates a design for maximum cavity impedance (shape parameter) 
even at a sacrifice in Qq. 


THE GAIN FUNCTION 

The gain function G^ (x) [1 + (x) ] summarizes the basic 

mechanism by which electrons orbiting through the cavity at the 
proper velocity, will drift into an accelerating or a retarding 
phase and hence, will absorb energy or give up energy to the 
cavity fields. 

A good physical description of the relativistic bunching 
mechanism, which dominates the interaction with a standing wave, 
has been given by Hsu (ref. 16). The key to this mechanism lies 
in the fact that any change in the kinetic energy of the orbiting 
electrons as they interact with the transverse electric field of 
the cavity, will primarily affect the radius of gyration. The 
gyrof requency is affected only slightly through the relativistic 
velocity dependence of the electron mass, and it is this slight 
change which can produce the desired bunching. 

The process can be illustrated by considering three elec- 
trons which enter the interaction region with the rotating elec- 
tric field vector in a decelerating (decel) , neutral, and accel- 
erating (accel) phase, respectively. The linearly polarized 
cavity field is here represented by two counter-rotating, circu- 
larly polarized vector components, only one of which will interact 
with the orbiting electrons. It is assumed that the gyrof requency 
of the entering electrons is set lower than the field frequency. 
Thus the first electron, which enters in the decel phase, would 
normally lag behind the rotating field vector and would eventu- 
ally fall into the accel phase. However, as the electron loses 
kinetic energy, this phase slippage is partially cancelled by 
the relativistic decrease in mass and the consequent increase in 
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gyrof requency , i.e., this electron tends to be "held" in the 
decel phase. 

On the other hand, the second electron which is initially 
in the neutral phase, will simply slip into the decel phase. 
Similarly, the third electron, initially in the accel phase, will 
gain in mass, decrease its gyrof requency , and will quickly also 
slip into the decel phase. The result is the desired bunching of 
electrons in the decel phase, i.e., a transfer of energy from the 
electrons to the cavity fields, or a net emission of microwave 
power. By similar reasoning, it is easy to show that if the 
gyrofrequency is set higher than the field frequency, net power 
absorption will result. Over a narrow range of gyrofrequency, 
then, the rotating electron cloud will go from a peak of positive 
to one of negative absorption, and this behavior has been veri- 
fied experimentally (ref. 2) . 

The gain function depends on the two variables: 


- 01 _ fl - 01 
k||U 0 " TT 


( 23 ) 


where ft= (e/m)B^ is the (relativistic) cyclotron frequency, and 


3 


k ll w 0 

oiu rt 


2 

TTWqT 


OIL 


( 24 ) 


A mono-energetic beam is assumed, with u = ug and w = wq . In 
effect, x describes the total phase lead or lag accumulated by 
an electron in passing through the cavity, referred to the phase 
of the rf electric field vector, whereas 3 is a measure of the 
(transverse) beam energy and of the transit time. 

A plot of G u (x) [1 + 3Q W (x) ] , with 


G w (x) 


2/ irx\ 
cos {—) 

(1 - x 2 ) 2 


( 25 ) 
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and 


Qoj(x) = x - 


, dG 

1 O) 

2G W dx 


0) 
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3c n c 


- x 


(26) 


is given in Figure 8, as calculated for the rectangular TEioi 
cavity (ref. 7). For each value of 6, there is a value of x 
(generally, x - -0.85) that will result in maximum negative 
absorption (beam loading) , and the amount of negative absorption 
increases as 6 is increased. Clearly, x - -0.85 is easily 
achieved by adjusting the magnetic field strength, but once more 
the maximum value of 6 is limited, for a given beam energy, by 
the physically realizable value of the transit time. A sharply 
peaked velocity distribution, therefore, emerges forcefully as 
an important design goal. 



Figure 8.- Gain function versus x, for the device 
shown in Figure 1. The parameter is 8. 
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The inverse dependence of Istart on the square of the 
transit time in effect defines an ideal condition where there is 
just a rotating electron cloud inside the cavity, with zero axial 
motion and, hence, infinite transit time. This condition cannot 
be reached in practice, of course, because of space-charge 
limitations, and because of the inevitable velocity spread on 
the beam, i.e., the slowest electrons will be turned around by 
the magnetic ramp (ref. 3), the fastest electrons will be out of 
synchronism, and only the electrons corresponding to the central 
part of the velocity distribution will actually be available for 
interaction with the cavity fields. 

A precise experimental check of the theoretical prediction 
of Istart i- s quite impossible until something about the actual 
velocity distribution of the beam and, hence, an average transit 
time, is known. One way to measure the velocity distribution is 
by means of a retarding-f ield analyzer. However, special care 
must be taken to allow for the non-rectilinear trajectory of the 
beam electrons (ref. 23) . 

A rough estimate of the transit time for electrons near the 
peak of the velocity distribution can be had by resorting to the 
gain characteristic. Clearly, oscillations generally will occur 
when the negative beam loading is maximum, so the slip parameter 
will have a value x = (ft -w/kjju) = -0.85, and with = (e/m)B, 
to, and k|| all known, this determines uq and, hence, the transit 
time. For example, with a 15-keV beam orbiting through a TE^o^ 
cavity resonant at 9.6 GHz, the axial energy for peak negative 
beam loading is 84 V, corresponding to a transit time of 4 nsec 
through the 21 mm cavity length. The kinetic energy thus is 
almost entirely transverse, and the relativistic mass change of 
the electrons is 2.93%. The value of Bi measured was 3439+10 
Gauss, and this corresponds to a relativistic cyclotron frequency 
that is 1.2 percent below the cavity oscillation frequency of 
9.463 GHz. 


One may go on to calculate the start-oscillation current as 
follows. The parameter 


8 


w. 


U) u. 


2.29 


and with this one enters the graph of gain function vs. slip 
parameter (Figure 8); this gives G w (1 + BQ U ) = -1.6. The shape 
parameter E/ioU = 3.06*10 6 ohms/m 2 for this cavity. Finally, 
from Eq. (15) one gets 


Istart^L = 75 • 3 roA (calculated). 



In comparing this value to the experimental result, one can, 
at best, bracket the "true" beam current within a maximum Cthe 
measured cathode current) and a minimum (the measured collector 
current). From reference 7, where the cathode current is used. 


I start Q L = 515 ^ ( measured ) / 


whereas, if the collector current were used, 

I , Q = 129 mA (measured) . 
start -L 


In either case, the agreement is better than order-of-magnitude . 
The main source of error in the theory probably is in the assump- 
tion of a delta-function velocity distribution (zero velocity 
spread) . 


CALCULATION OF EFFICIENCY 

The analysis of the previous sections was based on an 
expression for start-oscillation current that itself was derived 
from a linearized theory, with the assumption of a mono-energetic 
beam. The start-oscillation current probably is a reasonably 
good qualitative predictor of efficiency, in that it gives an 
indication of the intensity of the beam-cavity interaction. How- 
ever, as with any oscillator, the steady-state amplitude of the 
output is determined by non-linear, large signal phenomena, which 
call for entirely different computational procedures. Hsu 
(ref. 16) has formulated an analytical but very unwieldy expres- 
sion for the efficiency, based on including both first- and 
second-order terms in an assumed series solution to the Lorentz 
equation. However, in view of this early truncation of the series 
solution, it seems doubtful whether the resulting expression is 
really valid for large-signal conditions, and one is hardly 
tempted to go through the complex calculations necessary to derive 
a numerical result from Hsu's formulas. The Russian literature 
(refs. 4 and 6) gives estimates of the theoretical efficiency as 
high as 24 percent; however, as no derivations of their formulas 
are given, these estimates are difficult to accept without some 
skepticism. Thus, despite the obvious computational difficulties, 
a reliable estimate of the theoretical efficiency is probably best 
obtained by direct numerical integration of the electronic equa- 
tions of motion for a specific set of assumed conditions. 
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For completeness. Figure 9 summarizes the measured limiting- 
amplitude (saturated) behavior of the device shown in Figure 1, 
giving power output and efficiency in terms of the cathode cur- 
rent. The beam transmission to the collector was typically near 
25 percent, and a less conservative estimate of the efficiency, 
defined in terms of the minimum current known to be available 
for interaction with the cavity, would therefore be close to 
four percent at the maximum. 



Figure 9.- Power output and efficiency versus beam current 
for the device shown in Figure 1. 


A phenomenon which seems at first to be a serious limitation 
on the power output of cyclotron resonance oscillators, is the so- 
called cyclotron resonance breakdown (see Appendix B) . Fortun- 
ately, the margin between breakdown conditions and operating 
conditions can be made adequately wide, especially with high- 
velocity beams. 




CONCLUSIONS 


In the foregoing, results derived from the small-signal 
kinetic theory, together with some experimental data, have been 
used to bring out those physical parameters that dominate the 
electronic gain of cyclotron resonance devices. These are, the 
cavity geometry, the loaded Q, the transit time, and the single- 
electron gain function. The latter, in turn, is a function of 
the transverse kinetic energy, the synchronism conditions, and 
the transit time. A very long transit time is desirable, but is 
physically realizable only with a sharply peaked electron velocity 
distribution, and to achieve the latter, therefore, is an impor- 
tant design goal. A numerical example is given to indicate that 
the start-oscillation current can be predicted within a factor 
of 2-5 even with strong simplifying assumptions. 


Electronics Research Center 

National Aeronautics and Space Administration 
Cambridge, Massachusetts, January 1969 
109-02-03-01 
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APPENDIX A 


SPONTANEOUS RADIATION 


A well-known prediction of classical electromagnetic theory 
is that an accelerated charge will lose some of its kinetic 
energy by radiation (ref. Al) . For a single electron orbiting 
with radian frequency ft and tangential velocity w = 3c, the 
radiated power is 



If a large number of electrons are uniformly distributed in a 
circle spinning about its axis, the net radiation, however, will 
be zero. A uniform ring current does not radiate; some bunching 
around the circumference is necessary for net radiation to take 
place. If an ideal ensemble of identically phased orbiting 
electrons is assumed, one may linearly superimpose the radiated 
power contributions, and one may be tempted to use this as an 
estimate of the maximum conversion efficiency of a cyclotron 
resonance oscillator. It must be remembered, however, that such 
an estimate considers only_^the spontaneous cyclotron radiation 
and ignores the important J-E interaction with the cavity fields. 
Moreover, it is assumed that no energy is lost by radiation while 
the rotating electron cloud is being formed. Thus, the resulting 
expression probably has little validity in this case. 
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APPENDIX B 


CYCLOTRON RESONANCE BREAKDOWN 


Cyclotron resonance breakdown at low pressure in an evacua- 
ted microwave cavity , as discussed recently by Schrader (ref. Bl) 
has also been observed during our experiments (ref. B2) . Typi- 
cally, the cavity spontaneously changed its impedance to a 
milliwatt-level probing signal at the cavity resonant frequency 
as the ambient magnetic field was varied through cyclotron 
resonance. At pressures between 10“^ and 10 - 6 torr, breakdown 
typically occurred after a few microseconds. These results are 
in reasonable agreement with measurements by Schrader who found 
that the minimum power required for breakdown remains essentially 
constant at less than a milliwatt once the pressure is reduced 
below a certain minimum (in his case, 10“3 torr). He also found 
that the optimum magnetic field for breakdown is about 0.1 to 
0.2 percent above the value corresponding to rest-mass cyclotron 
resonance, as compared to 0.2 percent to 0.6 percent observed by 
the author. Physically, cyclotron resonance breakdown apparently 
is caused by trapping of electrons in the rf fields within the 
cavity, so that locally the effective pressure is increased to 
the point where ordinary high-frequency breakdown would be 
expected. Evidently, a breakdown of this type could seriously 
interfere with the operation of a cyclotron resonance oscillator. 
Fortunately, with reasonably energetic electrons, the magnetic 
field required to produce the appropriate relativistic cyclotron 
frequency 


n 


eB 

ym 0 ' 


where 
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1 



- 1/2 


r 


is 1 to 2 percent above the value that would produce breakdown of 
the evacuated cavity. A semi-empirical relationship is 


B ( start-osc , ) „ po 
B(cyclo. br.) ' 

and since this ratio should be made large, this is an argument in 
favor of high-velocity beams. 
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It is interesting to note that Antakov, et al. (ref. B3) , 
also report cyclotron resonance breakdown in a traveling-wave 
device using a cycloidal sheet beam. 
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